The Supernova Relic Neutrino Background 
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An upper bound to the supernova relic neutrino background from all past Type II supernovae is 
obtained using observations of the Universal metal enrichment history. We show that an unambigu- 
ous detection of these relic neutrinos by the Super-Kamiokande detector is unlikely. We also analyze 
the event rate in the Sudbury Neutrino Observatory (where coincident neutrons from v c D — > nne + 
might enhance background rejection), and arrive at the same conclusion. If the relic neutrino flux 
should be observed to exceed our upper bound and if the observations of the metal enrichment his- 
tory (for z < 1) are not in considerable error, then either the Type II supernova rate does not track 
the metal enrichment history or some mechanism may be responsible for transforming v^ T — * v e . 



I. INTRODUCTION 



A Type II supernova (SN n ) - the explosion triggered by the gravitational collapse of a single massive star - emits 
99% of its energy in neutrinos. The relic v c background created by all past SNu is potentially detectable in present 
and/or future large underground neutrino detectors. One of the goals of the current SuperK and SNO detectors is 
to detect this SRN background HIB] . The predicted SN n relic v c (SRN) flux depends crucially on the SN n rate as a 



to detect this SRN background |l|gl . The predicted SN n relic D c (SRN) flux depends crucially on the SN n rate as a 
\ function of redshift and the epoch of maximum SN rate (throughout the SN rate refers to the comoving densities of the 
J> ■ SNu rate) is important in determining the detectability of the SRN background. For example, if the SN rate should 
t— \ \ peak at redshifts of order 2 to 3, the majority of the SRNs would be rcdshiftcd to energies below typical detector 
0^ • threshold energies (~ 5 MeV). Since the same objects which are responsible for creating the SRN background are also 
responsible for the bulk of the heavy element production, knowledge of the SNu metal and neutrino production, in 
concert with the observationally inferred metal enrichment history of the Universe, is the straightest path (i.e., least 
model dependent) to predicting the flux and spectrum of relic v c from all past SN n []. It is our goal here to follow 
, this path in providing a generous upper bound to the expected SRN background. Furthermore, we account for the 
characteristics of the SuperK and SNO detectors, and use our calculated (upper bound to the) SRN background to 
predict (upper bounds to) the event rates at these detectors. These event rates are compared to expected backgrounds 
and to current limits. The current upper limit || (from the Kamiokande II detector) on the flux of supernova relic 
P e in the energy interval from 19 to 35 MeV is 226 cm _2 sec _1 . SNO is just beginning operation and has not decided 
+3 , upon a neutron detection strategy which is vital to the detection of v c . 

In the last few years progress has been made in constraining the recent star formation history of the Universe. In 
particular, a variety of observational evidence seem consistent with a comoving star formation rate (SFR) density 
which was much higher at redshifts z ~ 1 than at the present epoch. The history of star formation beyond z ~ 1 is 
less certain and it is not yet clear if the Universal SFR declined rapidly or evolved only mildly at higher redshifts. 
Support for this scenario comes from Pei & Fall || who have used chemical evolution models to explore the SFR and 
metal enrichment history inferred from observations of damped Lya systems. They find that the observed H i column 
densities may not represent the true column densities because of significant corrections due to dust. Since the Lya 
systems are identified from the spectra of quasars and since the Lya systems may contain dust, the implication is that 
some of the quasars may be invisible. This, in turn, suggests that some Lya systems may go undetected. When Pei 
& Fall correct for the effects of dust obscuration, they find evidence for rapid star formation at low redshifts (z ~ 1). 
This is to be compared to the predicted peak star formation rate epoch at redshifts of order 3-4 when obscuration 
is not taken into account. In particular, for their model with infall Pei & Fall find that the observational data is 
consistent with a SFR which increases until z ~ 1 and then decreases with further increases in redshift. Independent 
observations by Madau et al. [H of the metallicity enrichment rate (MER) are in excellent agreement with the Pei & 
Fall results. The direct quantitative support for Pei & Fall's model comes from the Canada- France redshift survey of 
faint galaxies (ll) which found that the comoving UV luminosity density of the Universe shows a sharp decline from 
z ~ 1 to the present. 



1 Type I supernovae are not expected to contribute appreciably to the relic neutrino background. 
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Here, using the model of Pei and Fall, we parameterize the metal enrichment history (MER) observed by Madau et 
al. jl4| to predict the SNu relic v e flux at Earth. In the past, similar calculations of this relic v e flux have been done 
by Totani & Sato |||, Totani, Sato & Yoshii Q, Malaney [|§ and Hartmann & Woosley Q. In contrast to almost 
all of the above analyses, we strive to minimize any model dependences by directly relating the supernova rate and 
its evolution to observations of the metal enrichment history to obtain supernova rate and, thereby, the SRN flux. In 
our calculations we make always make "conservative" choices of any uncertain parameters so as to obtain a robust 
upper bound to the SRN flux. From this upper bound we will conclude that it is unlikely SuperK will detect these 
relic neutrinos (because the signal will be buried under a large background event rate) and that the event rate in 
SNO should be vanishingly small. In §||, we outline the formalism for calculating the flux of SRN at Earth. In §111, 



we calculate conservative (i.e., generous) upper bounds to the relic v event rates at SuperK and SNO. In §[V 



wo 



review previous estimates of the SRN flux in comparison with ours and we examine neutrino oscillations as a possible 
mechanism for increasing the SN relic v c flux. 



II. THE SUPERNOVAE RELIC NEUTRINO SPECTRUM 



The spectrum of neutrinos at Earth due to all past supernovae depends on the differential (per unit energy interval) 
neutrino flux from each SN, on the redshift distribution of the SN rate, and on an assumed Friedmann-Robertson- 
Walker cosmology which may be parameterized by the Hubble parameter Hq and the matter density parameter f2o- 
For simplicity we ignore a possible cosmological constant at this point and discuss its effect later. If the supernova 
rate per unit comoving volume at redshift z is NgN(z) and the neutrino energy distribution at the source (at energy 
e) is Cfle), then the differential flux of relic neutrinos at Earth is given by 
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(1) 



where e' = (1 + z)e and the neutrinos are assumed to be massless. The angled brackets indicate that the dependence 
of the neutrino flux on supernova progenitors with different masses should be averaged over the initial mass function 
(IMF). In practice we will choose values for these average quantities so as to maximize the SRN background. 

The spectrum of the neutrinos from a supernova is parameterized as a Fermi-Dirac distribution with zero chem- 
ical potential, normalized to the total energy in a particular neutrino species (E„) emitted by the supernova, i.e., 
J Cf,(e)ede = E„. Then, for each neutrino species v, 
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The neutrino luminosity is thus characterized by E„ and T„ which, in turn, depend on the SN progenitor mass. 
However, the problem of obtaining the IMF-averaged neutrino flux simplifies because T„ docs not vary rapidly as the 
SN progenitor mass is changed . Adopting a flat Slo = 1 cosmology, and setting x = 1 + z , we can then write eq. |l| 
to a good approximation as 



(TV 



dx Ngpjfa;) 



exp(ea:/(T„)) + 1 



(3) 



where A = (120/7tt 4 ) cHq 1 = 1056 Mpc, with H = 50h 5Q km/s/Mpc. The results of Woosley et al. @ for a 
25 M Q supernova progenitor are used to fix (E„) = 11 x 10 52 ergs and (T„) = 5.3 MeV. The values of (E„) and 
(T„) characterize the detectable v e supernova neutrino spectrum . For comparison, recall that the data from SN 
1987A gave E u = 8 x 10 52 ergs and T„ = 4.8 MeV JJ]. Another issue concerns the role of the IMF and selection of 
the (Ej,) and (T„) values as averages. To obtain an upper bound to the detection rate, we use values of (E„) and 
(TV) which provide upper bounds to any reasonable average. This is possible because the flux integrated over the 
observable energy window (and hence the event rate) is an increasing function of both (E„) and (T„). The value 
of T„ is particularly insensitive to the progenitor mass because T„ derives its value from the temperature of the 
neutrinosphere formed during the collapse and the thermodynamic properties of the neutrinosphere (as long as it is 
well-defined) do not vary much with mass 0. In particular, T v ~ 4 — 5 MeV with 5.3 MeV being at the the upper 
end of the range. Thus, guided by models and SN 1987A, we parameterize the neutrino flux from supernovae so as 
to obtain a conservative upper bound to the SRN event rate. 

Under the assumption that the supernova rate tracks the metal enrichment rate, the supernova rate used to calculate 
the relic neutrino flux can be written as 
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where (Mz) is the average yield of "metals" (Z > 6) per supernova and pz is the metal enrichment rate per unit 
comoving volume. We have implicitly assumed that the metals come from SNn consistent with nucleosynthesis 
arguments || which show that the metal enrichment role of SN la is secondary to that of SNn . In any case, by 
neglecting SN la we overestimate the SRN flux (albeit, by only a factor of 2 at most). The other point to be noted 
here concerns our use of the metal enrichment history instead of the possibly more direct SFR to compute the SN 
rate. Both the SFR and the metal enrichment rate can be inferred from observations of the UV luminosity of star 
forming galaxies. But unlike the SFR which has a steep dependence on the adopted IMF, the metal enrichment rate 
is less sensitive to the IMF because the same (heavier) stars which are more UV luminous also eject more metals [[l4j . 
Thus, the SN rate is more closely tied to the metal enrichment history than to the SFR. 

To parameterize the evolution of pz(z) from the present back to z = 1, we use the results of Pei & Fall. In 
particular, we use the comoving metal production rate for their case with infall (Fig.l of Q) which is in good 
quantitative agreement with SFR observations at z < 1 Jn]-^|]. Since the neutrino flux from individual supernovae 
falls rapidly with increasing energy, and the lower energy neutrinos from high redshift supernovae are rcdshiftcd 
below the threshhold of detectability, our predictions are relatively insensitive to the high redshift (z > 1) behavior 
(as also noted by Hartmann & Woosley |3(J). This is fortunate since it is difficult to quantify precisely the z > 1 
evolution. For these reasons, we make the simplifying conservative assumption that the supernovae rate remains 
constant at higher redshifts: Nsn(z > 1) = Nsn(z = 1). It should be noted that the z < 1 evolution adopted 
here is likely quite robust in that independent studies reveal the same pattern of evolution (including, for example, 
that of the QSO luminosity density |l5[]) and the different observational data are in good quantitative agreement. 
Nevertheless, some changes to our adopted chemical enrichment history could be envisaged based on the arguments 
that the role of dust at high redshifts is still uncertain and, perhaps not all the star formation at higher redshifts has 
been observed |l6|-|l9|]. However, the relative insensitivity of our upper bound to the high redshift behavior insulates 
it against such uncertainty. 

To determine the average amount of metals ejected per supernova, the results of the calculations of supernova 
nucleosynthesis by Woosley and Weaver B are employed. From their published tables of the elemental composition 
of the ejecta, it can be ascertained that the heavy element yield ranges from Mz = 1.1 M for a 15 M Q SN progenitor 
to Mz = 4.2 M© for 25 M Q . These results are for an initial metallicity equal to 0.1 Z Q , which we assume characterizes 
the metallicity at redshifts around unity pi. In any case, the Mz values for SN progenitors with initial metallicity 
equal to Zq is greater by about 10-20% which, if used, would lead to a decrease in the predicted flux. Keeping in 
mind that the rate of events varies inversely as (Mz), we set (Mz) equal to 1 M© in the interest of obtaining an 
unambiguous upper bound. In Figure [l] we show the SRN spectrum that results from our adopted metal enrichment 
history and a conservative lower bound to the SN metallicity yields, (Mz) = 1 Mq. 



III. EVENT RATE AT SUPERK AND SNO 



It is not possible to detect SN relic neutrinos at all energies. For SuperK the observable energy window is likely to 
be from 19 to 35 MeV. Below 10 MeV, the 9 e from reactors [|o) and the Earth will completely overwhelm the relic 
neutrinos. Above 10 MeV and below the observable energy window, the main source of background is due to the solar 
neutrinos, radiation from outside the fiducial volume and spallation-produced events due to the cosmic-ray muons in 
the detector ||. Above 19 MeV the background is primarily due to atmospheric neutrinos |21jj . At energies greater 
than about 35 MeV, the rapidly (exponentially) falling SRN flux (peaked around 3 MeV) becomes smaller than the 
atmospheric D c flux, as can be verified from Figure |l|. Therefore the observable flux is obtained by integrating the 
differential flux over the neutrino energy range from 20.3 to 36.3 MeV (since e = E e + 1.3 MeV where E e is the 
energy of the positron and 1.3 MeV is the neutron - proton mass difference). We will also quote results in the more 
optimistic energy window of 15 - 35 MeV in the hope that with better background subtraction, SuperK will be able 
to probe these lower energy relic neutrinos. Detection of the SRN background in the much smaller SNO detector may 
be possible using coincident neutrons from v e D — > nne + . Because neutron detection at SNO is still in its infancy, we 
quote the total SNO event rate for positron energies above 10 MeV, corresponding to our SRN background. 

To calculate the event rate at SuperK, the detector is assumed to be 100% efficient in the observable energy window. 
The dominant reaction is v c p — > ne + with a cross section (<7 p (e)) two orders of magnitude larger than that of the 
scattering reaction [y Q e — ► v c e). The differential event rate in the interval de is then N p a p (e)j l/ (e)de and the predicted 
event rate at the detector is: 



R = -4N P pi r dx NarfaOVi C de t^-^t— ; , 
(T„) 4 J l J £1 exp(xe/(T v )) + l 
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where the q delineate the energy window (for this case, 20.3 and 36.3, respectively) and N p is the number of free 
protons in the detector. For SuperK, with a fiducial volume of 22.5 ktons, N p = 1.51 x 10 33 . 

Using the metal enrichment history to establish the supernova rate, the SN relic v e event rate at SuperK can be 
written as 

R = nnrfi (™®-\ f JM_\ (OA} events (R \ 

\(M Z ) J Vl0 52 er gS y \MeV J 22.5 kton- year 1 ' 

where we use a p (e) = 9.52 x 10~ 44 E e p e cm 2 [^2) with E e andp e (the energy and momentum of the positron) measured 
in MeV. We have set hso = 1 in the interest of obtaining an upper bound to the event rate. Also, for the same reason 
the average metal yield per supernova is taken to be 1 M Q , a lower bound to that obtained in the Woosley & Weaver 
models. For completeness we show in Figure [j] the differential rate of D c p — > ne + for our SRN background, with 
(E„) = 11 x 10 52 ergs and (T„) = 5.3 MeV. 
With our adopted SN parameters, the SRN event rate for a 22.5 kton-year exposure at SuperK is predicted to be 

R < 4 events 19 < E (MeV) < 35 , 

R < 7 events 15 < E c (MeV) < 35 . (7) 

Because the SRN spectrum falls rapidly with energy, the energy distribution of the events is strongly peaked at about 
10 MeV (see Fig. 1; in 5 MeV bins from 10 MeV to 40 MeV, the percentages are 37:29:17:10:5:2). If the threshold 
could be lowered to 10 MeV, our upper bound to the event rate at SuperK would increase to about 10/year. In terms 
of the flux at the detector, the results are as follows: the upper bound to the SRN flux integrated over all energies 
is 54 cm~ 2 sec _1 while in the relevant energy window from 19 to 35 MeV, the flux is 1.6 cm _2 sec _1 (to be compared 
to the current fjj upper bound of 226 cm _2 sec _1 ). In the larger energy window from 15 to 35 MeV, the observable 
flux is 3.7 cm~^3ec -1 . The reason for the large difference between the total and observable flux is two-fold. One, 
the observable energy window only captures the falling tail-end of the SRN spectrum and two, the event rate at low 
energies is artificially enhanced due to the SN rate which was assumed to be constant at high redshifts. 

In the energy window from 19 to 35 MeV, the expected background event rate from the atmospheric v c interacting 
with the protons (p e p — > ne + ) in the detector can be calculated. Using the atmospheric neutrino flux from Gaisser 
et al. the event rate for this background to the SRNs is only about 0.5/yr (for 22.5 ktons of water). However, 
there is another source of background which is dominant. The atmospheric muon neutrinos interacting with the 
nucleons (both free and bound) in the fiducial volume produce muons. If these muons are produced with energies 
below Cerenkov radiation threshold (kinetic energy less than 53 MeV), then they will not be detected, but their 
decay-produced electrons and positrons will. Consequently, the muon decay signal will mimic the v c p — > ne + process 
in SuperK. The event rate from these muon decays was estimated to be around unity for 0.58 kton-yr exposure of 
the Kamiokande II detector, forming the principal source of background after the various cuts had been implemented 
||. Extrapolating to the fiducial volume of 22.5 ktons for SuperK, we expect that SuperK should see ~39 events/yr 
as background to the SRN events. Although our predicted signal is much smaller than the sub-Cerenkov muon 
background, it may still be detectable because the energy distributions of the signal and the background are distinctly 
different. In such a case, a conservative criterion for the detectability of the signal is that it be greater than than the 
statistical fluctuations of the background. However, even with three years of data and assuming that the SRN flux is 
close to our upper bound, the SRN signal is only just about equal to the statistical fluctuations in the sub-Cerenkov 
muon background. This situation will improve, though not dramatically, if SuperK can lower its threshold (to SRN) 
to 15 MeV. 

Lastly we mention the SNO detector. Although much smaller than SuperK, the 1 kton SNO hopes to detect the 
SRN background by using the unique 2 neutron final state in v e D — > nne + . Using the cross section of Kubodera 
and Nozawa |p3[ , the upper bound to the event rate above 10 MeV is a not-very-promising 0.1/yr/kton. Again we 
show the differential event rate in Figure [l]. Note, however, that unlike the SRN signal in Super-K, this rate can be 
influenced by the large z SN rate, about which we know little. 



IV. DISCUSSION AND CONCLUSIONS 



A. Previous works 

Supernova relic neutrinos have been the focus of many previous studies [@,|0j2| -Q. The fluxes predicted in these 
studies spread over some two orders of magnitude, primarily due to the uncertain determinations of the present 
number density of galaxies, the SN rate in our galaxy at present, and/or the SN redshift distribution. More recently, 
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Totani et al. pg| used the population synthesis method to model the evolution of star-forming galaxies and they 
obtained a prediction for the flux of SRN. They found an event rate at SuperK (in the energy interval from 15 to 40 
MeV) of 1.2 yr _1 and the "most optimistic" prediction for their model was an event rate of 4.7/yr. Malaney used 
the Pei & Fall results in order to parameterize the evolution of the cosmic gas density which he then uses to calculate 
the star formation rate and, from that, the past supernova rate, finding a total SRN flux, integrated over all energies, 
of 2.0 — 5.4 cm~ 2 sec _1 depending on somewhat arbitrary low redshift corrections to the supernova rate. The work 
of Hartmann & Woosley |50|, using a SN rate proportional to (1 + z) 4 (motivated by J5 ll|)and normalized to the 
SN rate at present as derived from the Ha observations of the local Universe is most similar to ours. Their "best" 
estimate of a relic neutrino flux is ~ 0.2 cm _2 sec _1 . Although they do utilize Pei & Fall beyond z ~ 1, they (and 
we and others) have noted that this contribution is subdominant. However, Hartmann & Woosley f3C|| do not discuss 
the backgrounds to detecting the SRN and although their estimated flux is smaller than our upper bound by about a 
factor of five, they conclude the SRN may somehow be detectable. Although we agree with the Hartmann & Woosley 
estimate of the SRN flux in the sense that if we adopted their choices of parameters rather than our "conservative" 
choices we would predict the same flux, we disagree that this small flux is detectable. 

All these previous results are similar to, while less than the upper bound obtained in this paper. In fact, if we use 
our analysis of the SN rate alon g w ith the same IMF as used by Totani et al. for their spiral galaxies (which harbor 
most of the Type II supernovae |3l],^8|]), our SRN event rate in the 15 to 40 MeV range (for comparison with Totani 
et al. ) at SuperK falls to 1/yr. The total integrated flux falls to 11 cm _2 sec _1 while the result for the flux in the 15 
to 40 MeV energy window becomes 0.5 cm _2 sec _1 . These estimates agree well with those quoted from the previous 
works [Eq-pOl. In fact, the value of 1 event /yr obtained using the IMF from Totani et al. amounts to choosing the 
variables, (Mz), (T y ) and (E„) for an actual IMF rather than the extrema we have selected. It is more likely that 
any realistic IMF (chosen to fit other observables) when combined with a SN rate that peaks at z ~ 1 (as implied 
by the metal enrichment history) will yield an event rate that is an order of magnitude smaller than the upper bound 
we quote. Our upper-bound is robust because it is derived directly from the metal enrichment history which suggests 
that the SNu rate can peak no earlier than z ~ 1. 



B. Choice of Cosmology 



Throughout, we have assumed that f2 = 1 and q = 0.5. It is of some interest to ask how our SRN background 
predictions change if we change the background cosmology. Reducing the non-relativistic matter density from critical 
(fio < 1), by allowing positive curvature and/or a cosmological constant A, would reduce the expansion rate at late 
times and thereby increase the SRN flux, for the same Hq. The event rate increases by about 40% in going from an 
fio = 1 to an — 0.3, il\ — 0.7 Universe. But the estimation of the luminosity density (which is used to derive 
the metal enrichment rates) itself requires the assumption of a background cosmology and it typically increases less 
rapidly with redshift for cosmologies with smaller ilo [[ll| . These two effects tend to cancel out leaving the expected 
event rate nearly unchanged. Thus we do not expect our results to change substantially for a different background 
cosmology. 



C. Neutrino Oscillations 



The main goal of our work has been to obtain the most optimistic estimate of the SRN event rate at SuperK 
with the intent that if results from SuperK should exceed this upper bound, it could provide hints of new physics 
beyond the standard model. Here, we consider neutrino oscillations as a mechanism for maximizing the SRN flux. 
Since v x (where x= /i or r) only experience neutral current interactions, they decouple deeper in the SN where the 
temperature is higher. As a result, they stream out of the SN with a higher temperature than the v c . Because higher 
energy neutrinos are easier to detect, v c <-> D x oscillations have the potential to increase the SRN event rate. The 
maximum effect for any scenario is attained when the mixing is maximal. We will assume a mass hierarchy wherein 
the electron neutrino is the lightest (however, for an inverted mass hierarchy and resonant conversion in the presence 
of magnetic fields, see Q). This implies that the MSW resonance condition is not satisfied for v c <-» D x , but vacuum 
oscillations can still occur. If all three flavors are maximally mixed, then the oscillation probabilities average out to 
1/3 for any reasonable choice of mass differences because of the large distances traversed by the the relic neutrinos 
(typically of order of Hq • for the oscillation length to be comparable to this, Am 2 ~ 10 -25 eV 2 ). Such oscillations 
would make two-thirds of the original D c flux hotter as they would be "born" (would have oscillated from v x ) with the 
same temperature as the D x . To quantify the discussion here, we take (T Px ) = 2(T Pe ) (we might be exaggerating the 
spectral difference between P x and v c considerably here |33|] ) and assume that the same amount of energy is expelled 
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in all three flavors. This leads to an upper bound to the SRN event rate at SuperK of 11/yr with an observable flux 
of 4 cm _2 sec _1 in the 19 - 35 MeV energy window. The upper bound is larger by about a factor of 3 as a result 
of the increase in the number of neutrinos in the exponential tail of the neutrino distribution (where the observable 
energy window lies), due to the increase in temperature. As before, (Mz) has been set equal to 1 M©. For the case 
where v c is maximally mixed with only one of either or Sy, the upper bounds are 9/yr and 3 cm _2 sec _1 for the 
event rate and observable flux respectively. The upper bounds in the 15 - 35 MeV energy window are 14/yr for the 
three neutrino maximal mixing case, and 12/yr for the two neutrino maximal mixing case. In general, a decrease in 
the threshold (below 19 MeV) and neutrino oscillations seem to be required to boost the SRN flux to sufficient levels. 
Because the spectral shape of this oscillation enhanced SRN signal is sufficiently different from the sub-Cerenkov 
muon background, it may be detectable as a distortion of the expected muon background, if the SRN flux is in the 
vicinity of the upper bound we have quoted. For SNO, the two neutrino maximal mixing case gives an event rate of 
0.25/yr/kton while the three neutrino maximal mixing increases the rate to 0.29/yr/kton. 

A point to clarify here concerns the selection of the observable energy window given that the event rate in 19 - 35 
MeV window is now relatively large. Due to the larger signal, the relic neutrinos only become sub-dominant to the 
atmospheric neutrinos around 60 MeV (for a relic neutrino flux close to the upper bound). In fact, integrating out to 
a neutrino energy of order 60 MeV would increase the SRN event rate by about 50%. However the background from 
muon decay would increase by more than a factor of 3. One possibility this opens up is to also use the energy window 
from about 55 to 70 MeV since the muon-decay background is cut off at m M /2 (as decay occurs for muons at rest). 
For the same oscillation parameters as before, the upper bound to the event rate in the 55 to 70 MeV energy window 
is about 1/yr. However, the event rate due to the atmospheric neutrinos in the same energy window is of comparable 
magnitude and so, the situation is still not promising. 



D. Conclusions 



Using those observations most closely connected to the metal enrichment history of the Universe in order to relate 
the MER to the SNR, we have derived a robust upper bound to the supernova relic neutrino events at SuperK: 4 
events in the energy window from 19 to 35 MeV for a 22.5 kton-yr exposure. We have argued that the SuperK signal 
is dominated by SN n from z < 1 and so it is insensitive to the high redshift behavior of the metal enrichment rate. 
We use only the generic features of gravitational collapse SN models which have been substantiated by observations of 
SN 1987A to characterize the T> e spectrum emergent from SN n . In combination, these facts argue for the robustness 
of the upper bound to the SRN event rate obtained here. In addition, we have analyzed the backgrounds to the 
SRN events and conclude it is unlikely that SuperK will be able to detect these SRN neutrinos, unless the Type II 
supernova rate does not track the metal enrichment rate, or the observations of the star formation rate which lead to 
estimates of the metal enrichment rate at z < 1 are in considerable error, and/or some physics beyond the standard 
model is at play. We also find that the event rate at SNO will most likely be too small to be detected. The effect 
of flavor oscillations on the SRN flux has also been studied and the maximum possible increase in the event rate is 
less than a factor of 3. If the original flux is close to the upper bound quoted here and the mixing close to maximal, 
SuperK just might see the SRN flux as a distortion in its sub-Cerenkov muon background. 
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FIG. 1. The relevant event rates for the detection of supernova relic neutrinos (SRN), along with the predicted SRN spectral 
flux. The v e p rate is for SRN detection at SuperK, while the u e D rate applies for the SNO detector. The abscissa, energy, 
refers to the v a energy for all the cases except for the fx decay rate where it corresponds to the decay-produced electron's energy 
(E e ) plus 1.3 MeV (i.e., what the energy of a v c would have to be, in order to produce a positron with energy E c by 1> C + p 
reaction.) 
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